Examination of pencil beam predictions of small-field dose using dosimetry film and Monte Carlo simulations

The iPlan treatment planning sys-tem uses a pencil beam algorithm, with density cor-rections, to predict the doses delivered by very small (stereotactic) radiotherapy fields. This study tests the accuracy of dose predictions made by iPlan, for small-field treatments delivered to a planar solid wa-ter phantom and to heterogeneous human tissue using the BrainLAB m3 micro-multileaf collimator.

Relative biological effectiveness (RBE) and out-of-field cell survival responses to passive scattering and pencil beam scanning proton beam deliveries

The relative biological effectiveness (RBE) of passive scattered (PS) and pencil beam scanned (PBS) proton beam delivery techniques for uniform beam configurations was determined by clonogenic survival. The radiobiological impact of modulated beam configurations on cell survival occurring in- or out-of-field for both delivery techniques was determined with intercellular communication intact or physically inhibited. Cell survival responses were compared to those observed using a 6 MV photon beam produced with a linear accelerator. DU-145 cells showed no significant difference in survival response to proton beams delivered by PS and PBS or 6 MV photons taking into account a RBE of 1.1 for protons at the centre of the spread out Bragg peak. Significant out-of-field effects similar to those observed for 6 MV photons were observed for both PS and PBS proton deliveries with cell survival decreasing to 50-60% survival for scattered doses of 0.05 and 0.03 Gy for passive scattered and pencil beam scanned beams respectively. The observed out-of-field responses were shown to be dependent on intercellular communication between the in-and out-of-field cell populations. These data demonstrate, for the first time, a similar RBE between passive and actively scanned proton beams and confirm that out-of-field effects may be important determinants of cell survival following exposure to modulated photon and proton fields

Investigating the Implications of a Variable RBE on Proton Dose Fractionation Across a Clinical Pencil Beam Scanned Spread-Out Bragg Peak

Purpose: To investigate the clinical implications of a variable relative biological effectiveness (RBE) on proton dose fractionation. Using acute exposures, the current clinical adoption of a generic, constant cell killing RBE has been shown to underestimate the effect of the sharp increase in linear energy transfer (LET) in the distal regions of the spread-out Bragg peak (SOBP). However, experimental data for the impact of dose fractionation in such scenarios are still limited.

Methods and Materials: Human fibroblasts (AG01522) at 4 key depth positions on a clinical SOBP of maximum energy 219.65 MeV were subjected to various fractionation regimens with an interfraction period of 24 hours at Proton Therapy Center in Prague, Czech Republic. Cell killing RBE variations were measured using standard clonogenic assays and were further validated using Monte Carlo simulations and parameterized using a linear quadratic formalism.

Results: Significant variations in the cell killing RBE for fractionated exposures along the proton dose profile were observed. RBE increased sharply toward the distal position, corresponding to a reduction in cell sparing effectiveness of fractionated proton exposures at higher LET. The effect was more pronounced at smaller doses per fraction. Experimental survival fractions were adequately predicted using a linear quadratic formalism assuming full repair between fractions. Data were also used to validate a parameterized variable RBE model based on linear α parameter response with LET that showed considerable deviations from clinically predicted isoeffective fractionation regimens.

Conclusions: The RBE-weighted absorbed dose calculated using the clinically adopted generic RBE of 1.1 significantly underestimates the biological effective dose from variable RBE, particularly in fractionation regimens with low doses per fraction. Coupled with an increase in effective range in fractionated exposures, our study provides an RBE dataset that can be used by the modeling community for the optimization of fractionated proton therapy.

Differences between dual X-ray absorptiometry using pencil beam and fan beam modes and their determinants in vivo and in vitro

The aim of this study was to determine the influence of individual factors on differences in bone mineral density (BMD) using dual X-ray absorptiometry pencil beam (PB) and fan beam (FB) modes in vivo and in vitro. PB.BMD and FB.BMD of 63 normal Caucasian females ages 21-80 yr were measured at the lumbar spine and hip. Residuals of the FB/PB regression were used to assess the impact of height, weight, adiposity index (AI) (= weight/height(3/2)), back tissue thickness, and PB.BMD, respectively, on FB/PB difference. The Hologic Anthropomorphic Spine Phantom (ASP) was measured using the PB and FB modes at two different levels to assess the impact of scanning mode and focus distance. The European Spine Phantom (ESP) prototype, a geometrically well-defined phantom with known vertebral densities, was measured using PB and FB modes and analyzed manually to determine the impact of bone density on FB/PB difference and automatically to determine the impact of edge detection on FB/PB difference. Population BMD results were perfectly correlated, but significantly overestimated by 1.5% at the lumbar spine and underestimated by 0.7% at the neck, 1.8% at the trochanter, and 2.0% at the total hip, respectively, when using the FB compared with PB mode. At the lumbar spine, the FB/PB residual correlated negatively with height (r = 0.34, p < 0.01) and PB.BMD (r = 0.48, p <: 0. 0001) and positively with AI (r = 0.26, p < 0.05). At the hip, residual of trochanter correlated positively with weight (r = 0.36, p < 0.01) and AI (r = 0.36, p < 0.01). The FB mode significantly increased ASP BMD by 0.7% compared with PB. Using the FB mode, increasing focus distance significantly (p < 0.001) decreased area and bone mineral content, but not BMD. By contrast, increasing focus distance significantly decreased PB.BMD by 0.7%. With the ESP, the PB mode supplied accurate projected are of the bone (AREA) results but significant underestimation of specified BMD in the manual analysis. The FB mode significantly underestimated PB. AREA by 2.9% but fitted specified BMD quite well. FB/PB overestimation was larger for the low-density (+8.7%) than for the high-density vertebra (+4. 9%). The automated analysis resulted in more than 14% underestimation of PB. AREA (low-density vertebra) and an almost 13% overestimation of PB.BMD (high-density vertebra) using FB. In conclusion, FB and PB measurements are highly correlated at the lumbar spine and hip with small but significant BMD differences related to height, adiposity, and BMD. In clinical practice, it can be erroneous to switch from one method to another, especially in women with low bone density.

Long-term outcomes and prognostic factors of skull-base chondrosarcoma patients treated with pencil-beam scanning proton therapy at the Paul Scherrer Institute

BACKGROUND Skull-base chondrosarcoma (ChSa) is a rare disease, and the prognostication of this disease entity is ill defined. METHODS We assessed the long-term local control (LC) results, overall survival (OS), and prognostic factors of skull-base ChSa patients treated with pencil beam scanning proton therapy (PBS PT). Seventy-seven (male, 35; 46%) patients with histologically confirmed ChSa were treated at the Paul Scherrer Institute. Median age was 38.9 years (range, 10.2-70.0y). Median delivered dose was 70.0 GyRBE (range, 64.0-76.0 GyRBE). LC, OS, and toxicity-free survival (TFS) rates were calculated using the Kaplan Meier method. RESULTS After a mean follow-up of 69.2 months (range, 4.6-190.8 mo), 6 local (7.8%) failures were observed, 2 of which were late failures. Five (6.5%) patients died. The actuarial 8-year LC and OS were 89.7% and 93.5%, respectively. Tumor volume > 25 cm(3) (P = .02), brainstem/optic apparatus compression at the time of PT (P = .04) and age >30 years (P = .08) were associated with lower rates of LC. High-grade (≥3) radiation-induced toxicity was observed in 6 (7.8%) patients. The 8-year high-grade TFS was 90.8%. A higher rate of high-grade toxicity was observed for older patients (P = .073), those with larger tumor volume (P = .069), and those treated with 5 weekly fractions (P = .069). CONCLUSIONS This is the largest PT series reporting the outcome of patients with low-grade ChSa of the skull base treated with PBS only. Our data indicate that protons are both safe and effective. Tumor volume, brainstem/optic apparatus compression, and age were prognosticators of local failures.

Pencil Beam Scanning Proton Therapy for Pediatric Parameningeal Rhabdomyosarcomas: Clinical Outcome of Patients Treated at the Paul Scherrer Institute

BACKGROUND Parameningeal rhabdomyosarcomas (PM-RMSs) represent approximately 25% of all rhabdomyosarcoma (RMS) cases. These tumors are associated with early recurrence and poor prognosis. This study assessed the clinical outcome and late toxicity of pencil beam scanning (PBS) proton therapy (PT) in the treatment of children with PM-RMS. PROCEDURES Thirty-nine children with PM-RMS received neoadjuvant chemotherapy followed by PBS-PT at the Paul Scherrer Institute, with concomitant chemotherapy. The median age was 5.8 years (range, 1.2-16.1). Due to young age, 25 patients (64%) required general anesthesia during PT. The median time from the start of chemotherapy to PT was 13 weeks (range, 3-23 weeks). Median prescription dose was 54 Gy (relative biologic effectiveness, RBE). RESULTS With a mean follow-up of 41 months (range, 9-106 months), 10 patients failed. The actuarial 5-year progression-free survival (PFS) was 72% (95% CI, 67-94%) and the 5-year overall survival was 73% (95% CI, 69-96%). On univariate analysis, a delay in the initiation of PT (>13 weeks) was a significant detrimental factor for PFS. Three (8%) patients presented with grade 3 radiation-induced toxicity. The estimated actuarial 5-year toxicity ≥grade 3 free survival was 95% (95% CI, 94-96%). CONCLUSIONS Our data contribute to the growing body of evidence demonstrating the safety and effectiveness of PT for pediatric patients with PM-RMS. These preliminary results are encouraging and in line with other combined proton-photon and photons series; observed toxicity was acceptable.

Robustness of the Voluntary Breath-Hold Approach for the Treatment of Peripheral Lung Tumors Using Hypofractionated Pencil Beam Scanning Proton Therapy.

PURPOSE The safe clinical implementation of pencil beam scanning (PBS) proton therapy for lung tumors is complicated by the delivery uncertainties caused by breathing motion. The purpose of this feasibility study was to investigate whether a voluntary breath-hold technique could limit the delivery uncertainties resulting from interfractional motion. METHODS AND MATERIALS Data from 15 patients with peripheral lung tumors previously treated with stereotactic radiation therapy were included in this study. The patients had 1 computed tomographic (CT) scan in voluntary breath-hold acquired before treatment and 3 scans during the treatment course. PBS proton treatment plans with 2 fields (2F) and 3 fields (3F), respectively, were calculated based on the planning CT scan and subsequently recalculated on the 3 repeated CT scans. Recalculated plans were considered robust if the V95% (volume receiving ≥95% of the prescribed dose) of the gross target volume (GTV) was within 5% of what was expected from the planning CT data throughout the simulated treatment. RESULTS A total of 14/15 simulated treatments for both 2F and 3F met the robustness criteria. Reduced V95% was associated with baseline shifts (2F, P=.056; 3F, P=.008) and tumor size (2F, P=.025; 3F, P=.025). Smaller tumors with large baseline shifts were also at risk for reduced V95% (interaction term baseline/size: 2F, P=.005; 3F, P=.002). CONCLUSIONS The breath-hold approach is a realistic clinical option for treating lung tumors with PBS proton therapy. Potential risk factors for reduced V95% are small targets in combination with large baseline shifts. On the basis of these results, the baseline shift of the tumor should be monitored (eg, through image guided therapy), and appropriate measures should be taken accordingly. The intrafractional motion needs to be investigated to confirm that the breath-hold approach is robust.

Pencil -beam redefinition algorithm dose calculations for electron therapy treatment planning

The electron pencil-beam redefinition algorithm (PBRA) of Shiu and Hogstrom has been developed for use in radiotherapy treatment planning (RTP). Earlier studies of Boyd and Hogstrom showed that the PBRA lacked an adequate incident beam model, that PBRA might require improved electron physics, and that no data existed which allowed adequate assessment of the PBRA-calculated dose accuracy in a heterogeneous medium such as one presented by patient anatomy. The hypothesis of this research was that by addressing the above issues the PBRA-calculated dose would be accurate to within 4% or 2 mm in regions of high dose gradients. A secondary electron source was added to the PBRA to account for collimation-scattered electrons in the incident beam. Parameters of the dual-source model were determined from a minimal data set to allow ease of beam commissioning. Comparisons with measured data showed 3% or better dose accuracy in water within the field for cases where 4% accuracy was not previously achievable. A measured data set was developed that allowed an evaluation of PBRA in regions distal to localized heterogeneities. Geometries in the data set included irregular surfaces and high- and low-density internal heterogeneities. The data was estimated to have 1% precision and 2% agreement with accurate, benchmarked Monte Carlo (MC) code. PBRA electron transport was enhanced by modeling local pencil beam divergence. This required fundamental changes to the mathematics of electron transport (divPBRA). Evaluation of divPBRA with the measured data set showed marginal improvement in dose accuracy when compared to PBRA; however, 4% or 2mm accuracy was not achieved by either PBRA version for all data points. Finally, PBRA was evaluated clinically by comparing PBRA- and MC-calculated dose distributions using site-specific patient RTP data. Results show PBRA did not agree with MC to within 4% or 2mm in a small fraction (<3%) of the irradiated volume. Although the hypothesis of the research was shown to be false, the minor dose inaccuracies should have little or no impact on RTP decisions or patient outcome. Therefore, given ease of beam commissioning, documentation of accuracy, and calculational speed, the PBRA should be considered a practical tool for clinical use. ^

Characterization of the dose distribution in the halo region of a clinical proton pencil beam using emulsion film detectors

Proton therapy is a high precision technique in cancer radiation therapy which allows irradiating the tumor with minimal damage to the surrounding healthy tissues. Pencil beam scanning is the most advanced dose distribution technique and it is based on a variable energy beam of a few millimeters FWHM which is moved to cover the target volume. Due to spurious effects of the accelerator, of dose distribution system and to the unavoidable scattering inside the patient's body, the pencil beam is surrounded by a halo that produces a peripheral dose. To assess this issue, nuclear emulsion films interleaved with tissue equivalent material were used for the first time to characterize the beam in the halo region and to experimentally evaluate the corresponding dose. The high-precision tracking performance of the emulsion films allowed studying the angular distribution of the protons in the halo. Measurements with this technique were performed on the clinical beam of the Gantry1 at the Paul Scherrer Institute. Proton tracks were identified in the emulsion films and the track density was studied at several depths. The corresponding dose was assessed by Monte Carlo simulations and the dose profile was obtained as a function of the distance from the center of the beam spot.

Three-dimensional electron beam dose calculations

The MDAH pencil-beam algorithm developed by Hogstrom et al (1981) has been widely used in clinics for electron beam dose calculations for radiotherapy treatment planning. The primary objective of this research was to address several deficiencies of that algorithm and to develop an enhanced version. Two enhancements have been incorporated into the pencil-beam algorithm; one models fluence rather than planar fluence, and the other models the bremsstrahlung dose using measured beam data. Comparisons of the resulting calculated dose distributions with measured dose distributions for several test phantoms have been made. From these results it is concluded (1) that the fluence-based algorithm is more accurate to use for the dose calculation in an inhomogeneous slab phantom, and (2) the fluence-based calculation provides only a limited improvement to the accuracy the calculated dose in the region just downstream of the lateral edge of an inhomogeneity. The source of the latter inaccuracy is believed primarily due to assumptions made in the pencil beam's modeling of the complex phantom or patient geometry.^ A pencil-beam redefinition model was developed for the calculation of electron beam dose distributions in three dimensions. The primary aim of this redefinition model was to solve the dosimetry problem presented by deep inhomogeneities, which was the major deficiency of the enhanced version of the MDAH pencil-beam algorithm. The pencil-beam redefinition model is based on the theory of electron transport by redefining the pencil beams at each layer of the medium. The unique approach of this model is that all the physical parameters of a given pencil beam are characterized for multiple energy bins. Comparisons of the calculated dose distributions with measured dose distributions for a homogeneous water phantom and for phantoms with deep inhomogeneities have been made. From these results it is concluded that the redefinition algorithm is superior to the conventional, fluence-based, pencil-beam algorithm, especially in predicting the dose distribution downstream of a local inhomogeneity. The accuracy of this algorithm appears sufficient for clinical use, and the algorithm is structured for future expansion of the physical model if required for site specific treatment planning problems. ^

Improving 4D plan quality for PBS-based liver tumour treatments by combining online image guided beam gating with rescanning.

Pencil beam scanned (PBS) proton therapy has many advantages over conventional radiotherapy, but its effectiveness for treating mobile tumours remains questionable. Gating dose delivery to the breathing pattern is a well-developed method in conventional radiotherapy for mitigating tumour-motion, but its clinical efficiency for PBS proton therapy is not yet well documented. In this study, the dosimetric benefits and the treatment efficiency of beam gating for PBS proton therapy has been comprehensively evaluated. A series of dedicated 4D dose calculations (4DDC) have been performed on 9 different 4DCT(MRI) liver data sets, which give realistic 4DCT extracting motion information from 4DMRI. The value of 4DCT(MRI) is its capability of providing not only patient geometries and deformable breathing characteristics, but also includes variations in the breathing patterns between breathing cycles. In order to monitor target motion and derive a gating signal, we simulate time-resolved beams' eye view (BEV) x-ray images as an online motion surrogate. 4DDCs have been performed using three amplitude-based gating window sizes (10/5/3 mm) with motion surrogates derived from either pre-implanted fiducial markers or the diaphragm. In addition, gating has also been simulated in combination with up to 19 times rescanning using either volumetric or layered approaches. The quality of the resulting 4DDC plans has been quantified in terms of the plan homogeneity index (HI), total treatment time and duty cycle. Results show that neither beam gating nor rescanning alone can fully retrieve the plan homogeneity of the static reference plan. Especially for variable breathing patterns, reductions of the effective duty cycle to as low as 10% have been observed with the smallest gating rescanning window (3 mm), implying that gating on its own for such cases would result in much longer treatment times. In addition, when rescanning is applied on its own, large differences between volumetric and layered rescanning have been observed as a function of increasing number of re-scans. However, once gating and rescanning is combined, HI to within 2% of the static plan could be achieved in the clinical target volume, with only moderately prolonged treatment times, irrespective of the rescanning strategy used. Moreover, these results are independent of the motion surrogate used. In conclusion, our results suggest image guided beam gating, combined with rescanning, is a feasible, effective and efficient motion mitigation approach for PBS-based liver tumour treatments.

Segmented ionization chambers for beam monitoring in hadrontherapy.

Segmented ionization chambers represent a good solution to monitor the position, the intensity and the shape of ion beams in hadrontherapy. Pixel and strip chambers have been developed for both passive scattering and active scanning dose delivery systems. In particular, strip chambers are optimal for pencil beam scanning, allowing for spatial and time resolutions below 0.1 mm and 1 ms, respectively. The MATRIX pixel and the Strip Accurate Monitor for Beam Applications (SAMBA) detectors are described in this paper together with the results of several beam tests and industrial developments based on these prototypes.

An investigation of dose calculation accuracy in intensity-modulated radiotherapy of sites in the head & neck

Knowledge of the accuracy of dose calculations in intensity-modulated radiotherapy of the head and neck is essential for clinical confidence in these highly conformal treatments. High dose gradients are frequently placed very close to critical structures, such as the spinal cord, and good coverage of complex shaped nodal target volumes is important for long term-local control. A phantom study is presented comparing the performance of standard clinical pencil-beam and collapsed-cone dose algorithms to Monte Carlo calculation and three-dimensional gel dosimetry measurement. All calculations and measurements are normalized to the median dose in the primary planning target volume, making this a purely relative study. The phantom simulates tissue, air and bone for a typical neck section and is treated using an inverse-planned 5-field IMRT treatment, similar in character to clinically used class solutions. Results indicate that the pencil-beam algorithm fails to correctly model the relative dose distribution surrounding the air cavity, leading to an overestimate of the target coverage. The collapsed-cone and Monte Carlo results are very similar, indicating that the clinical collapsed-cone algorithm is perfectly sufficient for routine clinical use. The gel measurement shows generally good agreement with the collapsed-cone and Monte Carlo calculated dose, particularly in the spinal cord dose and nodal target coverage, thus giving greater confidence in the use of this class solution.

Investigation of stereotactic radiotherapy dose using dosimetry film and Monte Carlo simulations

This study uses dosimetry film measurements and Monte Carlo simulations to investigate the accuracy of type-a (pencil-beam) dose calculations for predicting the radiation doses delivered during stereotactic radiotherapy treatments of the brain. It is shown that when evaluating doses in a water phantom, the type-a algorithm provides dose predictions which are accurate to within clinically relevant criteria, gamma(3%,3mm), but these predictions are nonetheless subtly different from the results of evaluating doses from the same fields using radiochromic film and Monte Carlo simulations. An analysis of a clinical meningioma treatment suggests that when predicting stereotactic radiotherapy doses to the brain, the inaccuracies of the type-a algorithm can be exacerbated by inadequate evaluation of the effects of nearby bone or air, resulting in dose differences of up to 10% for individual fields. The results of this study indicate the possible advantage of using Monte Carlo calculations, as well as measurements with high-spatial resolution media, to verify type-a predictions of dose delivered in cranial treatments.